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We study gravitational collapse of the general spherically symmetric null strange quark fluid 
having the equation of state, p = (p — 4 B)/n, where B is the bag constant. An interesting feature 
that emerges is that the initial data set giving rise to naked singularity in the Vaidya collapse of 
null fluid gets covered due to the presence of strange quark matter component. Its implication to 
the Cosmic Censorship Conjecture is discussed. 
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I. INTRODUCTION 

A gravitating object when it undergoes indefinite collapse, the end product is singularity which is marked by 
divergence of physical parameters like energy density. As the singularity is approached, density diverges and it would 
therefore be of relevance to consider the state of matter at ultra high density beyond the nuclear matter. One of such 
possible states could be the strange quark matter which consists of u , d and s quarks. It is energetically most favored 
state of baryon matter. It could either be produced in quark-hadron phase transition in the early Universe or at ultra 
high energy neutron stars converting into strange (matter) stars JIJ. In the context of gravitational collapse, which is 
our concern here, it is the latter process which would be pertinent. 

The key question in the collapse process is whether the singularity so formed will be visible or will be covered by 
an event horizon prohibiting its visibility to an external observer. That the latter is the case goes by the name of 
Cosmic Censorship Conjecture (CCC) (see | 2 ] for reviews of the CCC)), which remains one of the most important 
unresolved issues in classical general relativity. In the strong form of CCC, there could emerge no null rays from 
the singularity in a reasonable space-time, and hence it is invisible for all observers. That is, there occurs no naked 
singularity for any observer. On the other hand the weak form states that null rays can emerge from the singularity 
which is however covered by an event horizon and hence they cannot reach out the external observer. In the weak 
form singularity is locally naked, say for an observer sitting on the collapsing star, but it is globally not because it is 
safely hidden behind an event horizon. 

There do however exist cases of regular initial data sets giving rise to possibility of null rays emanating from the 
singularity and reaching out to external observer (3). That is, both the possibilities of singularity being naked and 
hidden inside a black hole can occur. The critical question is what decides between these two possibilities? In this 
context, it has recently been argued that it is the shearing effect in the collapsing inhomogeneous dust cloud that 
is responsible for the ultimate outcome [Q. This happens because shear produces distortion in the collapsing fluid 
congruence which could cause distortion in the geometry of the apparent horizon surface. Such a distortion of the 
apparent horizon could let null rays emanating from the singularity to escape to external observer. It turns out that 
if shear close to the center exceeds a threshold limit, it gives rise to a naked singularity and else a black hole. 

Another critical question is, what is the state of matter as the singularity is approached? It is certainly a case of 
diverging density, and hence it would be appropriate to consider near the singularity matter in the highest known 
density form. That brings in the strange quark matter (SQM). Recently, collapse of charged strange quark fluid 
together with the Vaidya null radiation has been studied ||. In this paper, we would further like to analyze this 
process from the point of view of formation of naked singularity and its strength, and more importantly to bring out 
the effect of SQM. It turns out that effect of SQM leads to covering up the region in the initial data set window for 
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naked singularity. That is, it tends to favor black hole against naked singularity and consequently the CCC. This 
happens because SQM contributes an additional attractive potential. 

The SQM fluid is characterized by the equation of state p = {p — 45)/3 where B is the bag constant indicating 
the difference between the energy density of the perturbative and nonperturbative QCD vacuum, and p,p are the 
energy density and thermodynamic pressure of the quark matter |I], ||. The fluid consists of zero mass particles with 
the QCD corrections for trace anomaly and perturbative interactions [||. The boundary of a strange star is defined 
by P —> 0 which would imply p —> B. The typical value of the bag constant is of the order of B sa 10 15 g/cm 3 while 
the energy density, p ~ 5 x 10 15 g/cTO 3 0. This shows that SQM will always satisfy the energy conditions because 
p > p > 0. We shall however consider the equation of state p = (p — AB)/n as a generalization of SQM fluid, and 
particularly the cases n = 2, —> oo correspond to known cases of the Vaidya - de Sitter and the Vaidya in constant 
potential bath collapse respectively. 

In this paper, we shall first obtain the general solution for SQM null fluid with the generalized equation of state 
in the general spherically symmetric metric in the Bondi ingoing coordinates. We shall then bring out explicitly the 
effect of SQM on gravitational collapse in terms of covering of spectrum in the initial data set for naked singularity 
by finding the threshold values for the parameters involved in the mass function. 

The paper is organized as follows. In the next section, we obtain the general solution and analyze the collapse to 
show the effect of SQM in shrinking the parameter window in the initial data set giving rise to naked singularity. In 
the following section we discuss the strength of the singularity and we conclude with a discussion. 


II. STRANGE QUARK NULL FLUID COLLAPSE 


Expressed in terms of Eddington advanced time coordinate (ingoing coordinate) v, the metric of general spherically 
symmetric space-time [[ij 

ds 2 = —A(v, r) 2 /(u, r) dv 2 + 2 A(v, r ) dv dr + r 2 dfl 2 (1) 

dfl 2 = d0 2 + sin 2 9d(j) 2 . Here A is an arbitrary function. It is useful to introduce a local mass function m(v, r ) defined 
by / = 1 — 2 m(v,r)/r. For m = m(v) and / = 1, the metric (jl) reduces to the standard Vaidya metric. We wish 
to find the general solution of the Einstein equation for the matter field given by Eq.(|jj) for the metric (jl]), which 
contains two arbitrary functions. It is the field equation G 3 = 0 that leads to A(v,r) = g(v). This could be absorbed 
by writing dv = g{v)dv. Hence, without loss of generality, the metric (2) takes the form, 


ds 2 = - 


1 - 


2 m(v, r ) 


dv 2 + 2 dvdr + r 2 dVL 2 


The energy momentum tensor for the null fluid together with SQM can be written in the form 


( 2 ) 


Tab = pLk + (p + p)(l a n b + l b n a ) + pg a b 


( 3 ) 


Here p and p are functions of v and r and the two null vectors l a and n a 


— d a: Tl a — 


1 


2m(v, r) 


l a = (S“, n a = -5% - 


1 - 

1 f 2 m(y,r) 


A° _ A 1 

u a 

5'I 


Ll a = n a n a = 0 Ln a = -1. 


( 4 ) 


Substituting this in Eq. (||), we find Tg = T 3 = —p, T 2 = T| = p and Tg 1 = —p, and the trace T = —2 (p — p). Here 
p , p are the strange quark matter energy density and thermodynamic pressure while p is the energy density of the 
Vaidya null radiation. The Einstein field equations now take the form 











3 


At higher density, the equation of state becomes uncertain as is the case for nuclear matter and the strange quark 
matter would be no exception to it. It is therefore appropriate to keep the coefficient n free in the equation of state, 



P=-{p~ 4 B) 
n 

(6) 

for SQM. We should however be open to the possibility that in the unknown new super dense matter state, there 
could be altogether a different kind of contribution which could entirely change the situation. So far SQM is the most 
dense state of matter considered. 

Imposing the equation of state and combining Eqs. (|5b|) and (|5c|), we obtain the following differential equation 


u/ \ 2 , 327 tB 

m (v.r) = - m (v.r) H- r 

nr n 

(7) 

From this equation it is clear that the term involving the bag constant B makes the contribution similar to the 
cosmological constant. We shall thus seek the solution in the form, 


m(v,r) = mo(v, r) + —r 3 

O 

(8) 

which would lead to 

D _ (n + 1)A 

1671r 

(9) 

and the differential equation 

2 

mo(v,r) + — m' 0 (v,r) = 0 

(10) 

This has the general solution 

mo(v,r) = S(v)A n ~ 2)/n + M(y) 

(11) 


The two arbitrary functions M(v) and S(v) would be restricted only by the energy conditions. Here A is not the 
cosmological constant but instead is related to the bag constant B via Eq.([|). Thus the metric describing the radial 
collapse of null SQM in (y, r, 9 , <j>) coordinates reads as: 


ds 2 = —(1 — 


2 M(v) 2 S(v) 


Ar 2 , 

T 


- )dv 2 


2 dvdr + r 2 dfl 2 


( 12 ) 


^*2 j Tb 

This metric represents a solution to the Einstein equations for a collapsing null SQM. The physical quantities for this 
metric as in [|J || are given by 

1 


P = 


47 VI' 2 - 
1 


M{v) + S{v)r {n ~ 2)/n 


P = 


P = 


4tt r 2 

1 

4 n7rr 2 


n — 2 


n 


n — 2 


S(v) 


,-2/n 


A r 2 


S(v)r~ 2/n - A nr 2 


(13) 

(14) 

(15) 


Clearly, all the energy conditions would be satisfied for n > 2 because it would ensure p > 0 and p > 0, while p > 0 
would be taken care of when we choose the mass functions for both null radiation and SQM. The initial radius of the 
star from which the collapse begins would be given by p = 0 which would also relate the bag constant with the mass 
function S(v). 

We are studying the collapse of SQM null fluid on a non-flat but empty cavity. The first shell arrives at r = 0 
at time v = 0 and the final at v = T. A central singularity of growing mass would develop at r = 0. For v < 0, 
M(v) = S(v) = 0, i.e., we have 


A r 2 

ds 2 = — (1- 7 ^)dv 2 + 2 dvdr + r 2 dtt 2 


(16) 


and for v > T, M(v) = S(v) = 0, M = Mq > 0. The space-time for v = 0 to v = T is given by the generalized 
Vaidya metric and for v > T we have the generalized Schwarzschild metric: 


dg2= 2Mo _ Ar 2 2 

r 3 


2 dvdr + r 2 dtt 2 


(17) 
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A. Occurrence of naked singularities 

In this section, we adapt above solution to study the existence of a naked singularity. Let I\ a = dx a /dk be the 
tangent vector to the null geodesic, where k is an affine parameter. The geodesic equations, on using the null condition 
K a K a = 0, take the simple form 


d 2 v 1 \ M(y) 2 S(v) 

dk 2 r r r 2 / n 



(18) 


d 2 r 

dk? 


M(v) S{v ) 

ry* y* 2 jTl 



(19) 


Radial (6 and <fi = const.) null geodesics of the metric (11) must satisfy the null condition 

dr 1 
dv 2 

Clearly, the above differential equation has a singularity at r = 0, v = 0. If the singularity is naked, there must exist 
null ray emanating from it. By investigating the behavior of radial null geodesics near the singularity, it is therefore 
possible to determine whether outgoing null curves meet the singularity in the past. To go any further we would 
require specific form of functions M(y) and S(v), which we choose as follows: 

2 M(v) = av (a > 0), (21a) 

2 S(v) = f3v 2/n (fi > 0), (21b) 

Let X = v/r be the tangent to a possible outgoing geodesic from the singularity. In order to determine the nature 
of the limiting value of A' at r = 0, v = 0 on a singular geodesic, we let Xq = lim r ^o jj^o X = lim r ^o ?j->o -• Using 
©, © and L’HopitaPs rule we get 

oj 2 

Xq = lim X = lim — = lim — =-- (22) 

r ^0 v->o r >o v^o r r >0 v-o dr i_ a x 0 - (3X o /n - Ar 2 /3 


1 - 


2 M(v) 2 S(v 


A r 

r 2 / n 3 


( 20 ) 


which implies, 


aX 2 + /3Xo /n+1 ~X o + 2 = 0 (23) 

This is the equation which would ultimately decide the end state of collapse: a black hole or a naked singularity. 

Thus by analyzing this algebraic equation, the nature of the singularity can be determined. The central shell 
focusing singularity would atleast be locally naked (for brevity we have addressed it as naked throughout this paper), 
if Eq. (gH) admits one or more positive real roots 0. The values of the roots give the tangents of the escaping geodesics 
near the singularity. When there are no positive real roots to Eq. (^3|), there are no out going future directed null 
geodesics emanating from the singularity. Thus, the occurrence of positive roots would imply the violation of the 
strong CCC, though not necessarily of the weak form. Hence in the absence of positive real roots, the collapse will 
always lead to a black hole. The positive roots would define the range for the tangent slopes from which the null 
geodesics can escape to infinity. The critical slope would be given by the double root, marking the threshold between 
black hole and naked singularity. 


1. Case n = 3 

We now examine the condition for the occurrence of a naked singularity for n = 3. First note that the Eq. ( |23| ) 
is free of A and hence it has no effect on the question under study. However its presence makes the background 
space-time asymptotically non flat. This happens because when r —> 0 the term Ar 2 /3 in Eq.(|2^) tends to zero. For 
fi = 0, the allowed range for a is given by (0, 1/8] as obtained earlier ||To| for the Vaidya null radiation collapse. In 
this case, it would be black hole for a > 1/8. 

The numerical computation reveals that Eq. (^3|) would always admit two positive roots for a < ac■ Tangent to 
all outgoing radial null geodesics would lie in the range X 2 < A' < X±, where Ai and X 2 are the two roots. Table I 
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TABLE I: Variation of ac and Xq for various /3 (n = 3) 


(3 

Critical Valueac 

Equal Roots A'o 

0 

1/8 

3.9999 

0.05 

0.093728958525 

4.18083 

0.1 

0.06294108366 

4.39273 

0.15 

0.032689168719 

4.64701 

0.2 

0.0030431168445 

4.9625 


TABLE II: Values of Roots Xg for a < ac for different (3 (n = 3) 

13 

a < ac 

Roots (Ao) 

0.0 

0.1 

1.40338, 1.9342 

0.05 

0.09 

1.52319, 1.7237 

0.1 

0.06 

1.55259, 1.74615 

0.15 

0.03 

1.57943, 1.78409 

0.2 

0.003 

1.69165, 1.72026 


shows the critical values of a for various values of (3 . The window for naked singularity is defined by (0,ac], and it is 
black hole for a > ac- Table II indicates the slope range, given by the two roots, for the null geodesics to escape. It 
is seen that ac decreases with increase in /3, i.e., initial data set (0,1/8] for a naked singularity of the Vaidya collapse 
shrinks by the introduction of SQM. There exists a threshold value (3t = 0.205198 such that for (3 > (3t , gravitational 
collapse of strange quark null fluid would always end into a black hole for all a. 

Note that a refers to rate of collapse of the null radiation while (3 would refer to that of SQM. The /3-threshold 
would therefore define a critical rate of collapse for SQM required for collapse to end in a black hole. Then it would 
fully respect CCC. Thus introduction of quark matter favors formation of black hole. 


2. Other Cases 


a. Case n = 2 
Then we find that 


1 

P An rr 2 


M(v) + f3 , 


and the algebraic equation takes the form 



{a + /3)Xq — Aq + 2 — 0 


(24) 


(25) 


The metric in this case takes the form of the Vaidya - de Sitter metric. The singularity is visible for (a + (3) < 1/8. 
It is the null fluid collapse in the background of the de Sitter space, where A is generated by the bag constant. 

b. Case n —> oo. 

We have 


/X= W^ (W )’ P = 0, P = 


B 

Am ' 2 


(26) 


TABLE III: Variation of /3t with n 


2/n 

(3t 

0.9 

0.144 

0.8 

0.167 

0.75 

0.18 

0.5 

0.2052 

0.25 

0.45 

0.125 

0.6194 

0.0625 

0.756 
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TABLE IV: Values of equal roots A'o for different n and (3 


2/n 

0 

OLC 

Equal Root (Xo) 

0.90 

0.14 

0.00343411423125 

4.21544 

0.80 

0.16 

0.00502531901422 

4.47484 

0.75 

0.17 

0.00679411554354 

4.61836 

0.5 

0.20 

0.0306384999095 

5.1784 

0.25 

0.44 

0.00178952699837 

9.51405 

0.125 

0.61 

0.00077533872957 

16.519 

0.0625 

0.74 

0.00061702385153 

27.1684 


In this case we have the dual Vaidya metric or Vaidya metric with constant potential [|ll|. The algebraic equation: 
aXfi + (/3 — l)Xp + 2 = 0, would admit a positive root for a < l/8(/3 — l) 2 , giving the range for naked singularity as 
obtained in |ll| . This is simply the null fluid collapse in the background of constant potential which is characterized 
by Tq = Tj 1 = const./r 2 , as is the case in Eq. (p6|) above. Note that (3 < 1 else the metric signature would change, 
c. Other n 

We also note that as n increases, so does the threshold value /3 t- This is shown in Tables III and IV. 


III. CURVATURE STRENGTH OF SINGULARITIES 

An important aspect of a singularity is its gravitational strength |l^| . A singularity is gravitationally strong in the 
sense of Tipler, if it destroys any object which falls into it and weak otherwise. It is now widely believed that space- 
time does not admit an extension through a strong curvature singularity, i.e., space-time is geodesically incomplete. 
Through a weak singularity, space-time could be analytically extended to make it geodesically complete. There have 
been attempts to relate strength of a singularity to its stability |b| . Recently, Nolan [Q gave an alternative approach 
to check the nature of singularities without having to integrate the geodesics equations. It was shown |Q that a 
radial null geodesic which runs into r = 0 terminates in a gravitationally weak singularity if and only if dr/dk is finite 
in the limit as the singularity is approached (this occurs at k = 0, the over-dot here indicates differentiation along the 
geodesic). If the singularity is weak, we have 


dr 

dk 


d 0 k 


Using the asymptotic relationship ( dv/dk ~ doX 0 and v ~ doX 0 k ) and Eq. (pl|), the geodesic equation yields 

^- (aXodv 1 k~ 1 + fiX^d^k- 1 - ^d 0 k)d%X 2 


(27) 


(28) 


But this gives 


d 2 v 

dk? 


ck 1 


(29) 


where c = {aX^ 1 ^ 2 ^ n + P)X^ i ' n+1 ^ n d^ 1 . This is inconsistent with dv/dk ~ doX 0 , which is finite. Since the coefficient 
c is non-zero, the naked singularity is gravitationally strong in the sense of Tipler ]l5[ . Having seen that the naked 
singularity is a strong curvature singularity, we check it for scalar polynomial singularity. The Kretschmann scalar 
for the metric with the prescriptions (pl|), takes the form 


K = 


a 2 n 4 Xg + (4/3A n 2 - 6/3An 3 + 2/3An 4 )X 0 2/ " + (12a/?A n 2 + 


8a(3An 3 + 6a/3A n 4 )X^ n+2)/2 
+ (12 p 2 n + 15/3 2 n 2 + 3/3 2 n 4 )X 4 ' r 


-A 2 

3 


(30) 


which diverges at the naked singularity and hence the singularity is also a scalar polynomial singularity. The Ricci 
scalar also diverges. It however vanishes for the Vaidya space-time |[T7i| . Thus the naked singularities studied here are 
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strong curvature singularity and hence are physically significant. Lastly, we shall calculate Weyl scalar 


C = 


4 

3n 4 r 4 


an 1 2 3 4 X + 2/3X 2/n + 3 (3X 2/n + /3n 2 X 2/n 


(31) 


which too would diverge. The Weyl curvature describes non local effects of gravitation produced by free part of the 
field. It is generated by inhomogeneity and anisotropy, particularly divergence of shear Jl7| . In the context of naked 
singularity, like shear and inhomogeneity the Weyl curvature would also play significant role. 


IV. DISCUSSION 

In this paper, we have obtained the general solution for null SQM fluid with the equation of state given by Eq. (](]) 
for the general spherically symmetric metric (|^) in the Eddington advanced time coordinate (ingoing coordinates). 
We have used the solution to study the end state of the collapse. The present case is an example of non self-similarity 
as well as non asymptotic flatness, and yet there does occur a regular initial data set which would lead to naked 
singularity. 

The relevant question is what effect does the presence of the SQM have on formation or otherwise of a naked 
singularity. Our results imply that the presence of SQM leads to shrinking of the initial data space for naked 
singularity of the Vaidya null fluid collapse. That is, it tends to favor black hole. This tendency is caused by the 
additional attractive potential, varying as r -2 /™, produced by SQM which results in strengthening of gravity. There 
also exists a threshold value (3t = 0.205198 such that for (3 > (3t the end state of collapse of null SQM is always a 
black hole for all a. This is the critical value of the rate of collapse of SQM for respecting CCC. The case n = 3 
represents the zero mass particles with QCD corrections for trace anomaly and perturbative interactions Q. The 
energy conditions require n > 2. The case of n = 2 corresponds to the null fluid collapse in the background of de 
Sitter space where the bag constant provides the A. The case n —> oo, corresponds to the null fluid collapse in the 
background of constant potential space as studied in |ll| . These are the two extreme limiting cases encompassing the 
physically allowed cases. Though there is no much physical motivation in the context of SQM for n / 3 cases, yet 
these two particular cases are interesting. That is at least these three cases could be considered in a unified equation 
of state given by Eq. (8). 

As mentioned earlier the strong version of the CCC doesn’t allow even locally naked singularity, i.e., the space-time 
should be globally hyperbolic. It turns out that necessary and sufficient condition for a singularity to be locally naked 
is that the algebraic Eq. ( |23| ) should have atleast one or more positive root [|| . Hence existence of the positive roots 
of Eq. (H|) is a counter example to the strong version of the CCC. In the absence of the proof of any version of the 
CCC, such examples remain the only tool to study this important and unresolved problem. 

Quark stars could be formed in the realistic astrophysical setting. The core collapse of a massive star after the 
supernova explosion sets in first and second order phase transitions which result into deconfined quark matter. The 
other possibility is that some neutron stars could accrete matter and undergo phase transition to turn into quark stars 
fj~8| , p~9 . |20| . Thus study of gravitational collapse with quark matter component is quite in order because it is perhaps 
astrophysically more realistic. In the ultimate stage of collapse close to the singularity, density is diverging. The 
quark matter contribution would therefore perhaps be most significant in deciding the ultimate result of the collapse. 
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